Weak limits for quantum random walks 
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We formulate and prove a general weak limit theorem for quantum random walks in 
one and more dimensions. With X n denoting position at time re, we show that X n /n 
converges weakly as n — > oo to a certain distribution which is absolutely continuous 
and of bounded support. The proof is rigorous and makes use of Fourier transform 
methods. This approach simplifies and extends certain preceding derivations valid 
in one dimension that make use of combinatorial and path integral methods. 
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I. INTRODUCTION 

Let Ri, i?2, • • • be independent identically 
distributed random variables taking values in 
the reals R, and suppose that they have com- 
mon mean fi = E(R{) and finite non-zero 
variance a = E{R\) - /j, 2 . The central limit 
theorem asserts that the sum X n 
satisfies 



EILt Ri 



X n - n\x 



N as n — > oo 



where N denotes the normal (Gaussian) dis- 
tribution with mean and variance 1, and 
denotes weak convergence: 

T n ^T if E(f(T n ))^E(f(T)) (2) 

for all bounded continuous functions / : R — > 
R. An early version of this now classical the- 
orem for random walks was proved as long 
ago as 1733 by de Moivre, In the modern 
theory, the conditions on the Ri are relaxed 
to allow non-independent non-identically dis- 
tributed random variables taking values in 
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general spaces. Since the weak limit of X n , 
suitably normalized, depends only on the 
probability measures associated with the X n , 
we may think of the central limit theorem as 
a result about weak limits of measures, rather 
than about the stochastic process (X n : n > 
1) itself. This is an impoverishment of the 
theory, since it overlooks the random vari- 
ables themselves. 

There has been recent interest, [lj, in a 
new type of process termed a quantum ran- 
dom walk. Quantum random walks give rise 
to certain sequences (fi n : n > 1) of probabil- 
ity measures, each of which is given in terms 
of the preceding measures in the sequence. 
Whilst it is possible as always to construct 
random variables having these measures, this 
may not be done in a natural manner as in the 
theory of stochastic processes. One may nev- 
ertheless ask whether, subject to an appropri- 
ate normalization, the fi n converge weakly to 
some non-trivial distributional limit. Results 
in this direction have been obtained for one- 
dimensional quantum walks by Konno 0, HI . 
We show in this note how to simplify and 
extend such results. We introduce a new 
method of studying such weak limits, and we 
apply this method to quantum walks in one 
and higher dimensions. 
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We consider first a quantum random walk 
on the integers Z. At each time n (e N) the 
state of the particle is transformed by a uni- 
tary operator described by a rotation of the 
internal degree of freedom followed by a con- 
ditional shift of the position, pj; the internal 
degree of freedom represents a coin that de- 
termines the shift of the position. The over- 
all state of the system belongs to the Hilbert 
space H c <8> Hp, where He is associated with 
the internal degree of freedom (coin space) 
and Hp with position. In the simplest case, 
we have H c = C 2 and H P = £ 2 {Z). A suit- 
able basis for Hp is given by the eigenstates 
of the position operator X 



XV , 



x G Z, 



(3) 



subject to (v x ,v x >) = 5 XX >, the Kronecker 
delta. A general state of the system may be 
written with respect to this basis as 



<4) = ^2^2^ j (x)v x w j , 



(4) 



where the vectors Wj, j = 1, 2, define a stan- 
dard basis in H c . The probability fjL n (x) of 
finding the particle at the position x at time 
n is given by the standard rule 



(5) 



where ip n = U n ipo with U the time-evolution 
operator of the walk and ipo the initial state 
of the system. 

The asymptotic properties of the sequence 
{fJ>n '■ n > 1) are studied in the next sec- 
tion. Such results are extended in Section 
III to quantum walks in two and more di- 
mensions. We highlight two special features 
of such asymptotics, namely: instead of nor- 
malizing by \fn as in f^Q), we shall normalize 
by n, and the weak limit is absolutely contin- 
uous with bounded support. 



II. WEAK LIMIT FOR 
ONE-DIMENSIONAL QUANTUM 
WALKS 

In order to define the position of a quan- 
tum particle as a random variable, we con- 
sider the evolution of the position operator 
in the Heisenberg picture starting from time 
n = 0. At each time n, the eigenvalues of the 
operator X n = W n XU n define the possible 
values of the particle's position with corre- 
sponding probability given by (jSJ), where the 
dependence on the initial state ip is explicit. 

Although the position may be treated as 
an ordinary random variable, the sequence 
(X n : n > 1) does not define a stochastic 
process, since the simultaneous measurement 
of X n for different n would change the quan- 
tum random walk at each step. Therefore we 
let the system evolve repeatedly under U up 
to time n, without measuring it, and then we 
study the properties of the distribution jj, n of 
X n . 

Let ipo be any initial state in He® Hp with 
all moments E(X r ) finite. In order to sim- 
plify the calculations which follow, we con- 
sider transformations in terms of wave func- 
tion components, and we take the Fourier 
transform space £ 2 (Z) = L 2 (K), where K = 
[0, 2-71") is thought of as the unit circle in M 2 . 
We define an inner product on L 2 (K) by 



2n 



mm 



dk 
2^ 



(6) 



and we note the isometry between £ 2 (Z) and 
L 2 (K) given by 



(7) 



with inverse 



ixk- dk 



ip i— > ip where ip(x) = e ip(k) 

Jo 2tt 

(8) 

The right shift S on £ 2 (Z) given by 
S(i>x)™oc = (V^-O-oo corresponds to the 
multiplication operator Sip = e lk ip on L 2 (K). 
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Our fundamental Hilbert space is thus 
H = H c ® L 2 (K), the space of C 2 - valued 
functions 

on K satisfying 

W\* = \Mb + m\b < <*>■ 

As usual, we consider state vectors normal- 
ized by || ijj || 2 = 1. The evolution of the walk 
comprises repeated applications of an inter- 
nal transformation (coin toss) A acting on 
C 2 , followed by the shift S given by 

o(Mk)\_ ( e ik Mk) \ nm 
Thus the total evolution U on H is given by 

w -(C.-)H#))' Wk) ' 

(11) 



If we begin the quantum random walk with 
an initial state ^>q £ H, its state after n steps 
is 

V n = U n V = U(k) n V (k). (12) 

For each k, U(k) has two eigenvalues Xi(k) 
and \2(k) with |Aj(fc)| = 1, and has corre- 
sponding eigenvectors vi(k), V2(k) G C 2 that 
define a basis for H. We assume henceforth 
that 

Ai(fc) ^ A 2 (fc), (13) 

since otherwise U(k) is diagonal; by (fTTj) then 
A is diagonal and the state evolves trivially, 
either to the right or to the left. 

The mapping k h- > U(k) is C°° and the 
eigenvalues are distinct for each k, and there- 
fore the eigenvalues Xj(k) are C 00 functions of 
k, and the eigenvectors Vj(k) may be chosen 
to be C°° with normalization ||%-(A;)|| = 1. 
By expanding the wave function in terms of 
this basis, the nth time evolution becomes 



= U(k) n Mk) = Xi{k) n (vx(k), Mk))vi(k) + A 2 (fc)> 2 (fc), ^ {k))v 2 (k), (14) 



where each component on the right hand side 
is a C°° function of k. The moments of the 
position distribution are given in terms of the 
operator X according to the standard for- 
mula 

E(X r n ) = (^ n ,X r ^ n ). (15) 

Using the isometry between £ 2 (Z) and 
L 2 (K), the above expectation may be writ- 
ten as 

f 2n dk 
E(X:)= / {* n {k),D*9 n (k))£-, (16) 

where D = X = —id/dk is the position op- 



erator in the momentum space L 2 (K). For 
fixed r we can compute D r ^ n {k) by (fT3j) and 
Leibniz' rule. It is easily seen that 



j 

(17) 



where (n) r = n(n — 1) • • • (n — r + 1). Equa- 
tions (fTfijl and (|T7jl yield, as n — ► oo, 
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E[{X n /nY\ = / ^j(k) n ~ r (DXj(k)) r (vj(k), ty (k))(ty n (k),Vj(k)) - — h 0(n -1 ) 
Jo , 27r 



2?r 

E 



DXjjk) 



!(^(A;),vI/ (A;))| 2 ^ + O(n- 1 ). 



(18) 



Let f2 = Kx{l, 2}, let /i be the probability As an example, we consider some specific 

measure on Q given by | (^ (k), Vj(k)) \ 2 dk/2n cases of unitary quantum walks. We consider 

on K x {j}. Let hj(k) = Xj{k)~ 1 DXj{k) and first the Hadamard matrix 

define h : Q — > R by h(k,j) = hj(k). (h is i / \ 

real because |A 3 -(fc)l = 1.) By (HD, A = J_ f j _ j J . (21) 



£[pf n /n) r ] -> / /fd/i asn-^cx). (19) By simple calculus, 
in 



Since /i is bounded and the above relation 
holds for all integers r > 0, we deduce by the 
method of moments the following. (See [f| 
for the general theory of weak convergence.) 

Theorem 1. With notation as above, 

-X n Y = h{Z) asn^oo, (20) 
n 

where Z is a random element of Q with dis- 
tribution ji. 

In particular, the support of Y is 
[min h, max/t], the range of h, at least pro- 
vided the density of fi given above does not 
vanish on some interval. 

A similar weak limit theorem for X n /n 
has been proved by Konno 0, HI , by different 
methods and with a quite different descrip- 
tion of the limit. 

We note that no assumption has been 
made above on the matrix A and the initial 
state i/; , and thus the above result holds for 
any unitary quantum walk on the integers, 
subject only to (JTBJ). Note also that \i de- 
pends only on the overlap between the initial 
state of the system and the eigenvectors of 
U(k), whereas h depends only on the coin 
flip matrix A. 



Xj(k) = -j= sin k ± ^Jl — | siir /. 



and thus 



-iX' 3 {k) 



cos k 



2 - sin 2 k 



(22) 



(23) 



Hence the limit distribution is concentrated 
on the interval 



[min h, max h] 



1 1 

71' 71 



(24) 



For a general unbiased walk, we take as 
coin flip the unitary matrix 



I / e? 

= 71 -e"^^ J' (25) 

where (p,i/j G R, with corresponding evolu- 
tion 



U(k) 



e ik e -i((p-ip) 



U^{k). (26) 



With ip + ip = a, ip — ip = b, the eigenvalues 
may be written in the form 



Xj(k) = —j= sin(A; + a) ± y 1 — | sin 2 (A; + a), 

(27) 
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and therefore 



h{k,j) 



cos(/c + a) 



2 - sin (A; + a) 



(28) 

Thus the general unbiased walk has exactly 
the same behaviour as the Hadamard case, 
subject to a shift in the momentum parame- 
ter of the wave amplitudes. We have as before 
that the domain of the limit distribution is as 
in pijl. 

Finally we introduce a "biased" random 
walk by defining a bias factor p in the coin 
flip matrix 



U( P ) 



VP V 1 - p 



Vi - p ~Vp 

that gives rise to the evolution 



(29) 



U p (k) 



e ik sfp~ e ik y/l - p 
e~ ik y/l — p —e~ ik sfp 



(30) 



The evolution under U of a general two- 
component wave-function corresponds to 



■w 2 



e ik Q 

e~ ik 



( vpMk) + V^pMk) 



(31) 



where the two internal states transform dif- 
ferently. In fact, the first component receives 
a kick of momentum +k with probability p 
and —k with probability 1 — p; the opposite 
holds for the second component. 
In terms of p, the eigenvalues are 



Xj(k) = iy/psin k ± 
and thus 

-iX'^k) _ 



l-psin 2 k (32) 



h(k,j] 



cos k 



Hk) 



sin 2 k 



(33) 



on the asymptotic momentum distribution by 
changing the support of the limit distribu- 
tion. 

The representation (|2(Jj) of the limit vari- 
able allows a direct computation of the 
asymptotic probability density function in 
most cases of interest. For example, assume 
that the initial state is at position 0. If the 
coin initially is in a given state i — 1 or 2, 
then ^o(k) = (q) or (5), respectively, and 
thus p = \vji(k)\ 2 dk/2ir on K x {j}. If we 
consider instead a random initial state of the 
coin, we have a mixture of these two pure 
states and thus 



1 2 

-J2\ v d k )\ 



,dk 
27 



dk 
Att 



(34) 



on IK x {j}; that is, p is the uniform dis- 
tribution on Q. In the Hadamard case, for 
example, with h given by (f2~3*j) . if X = and 
the coin initially random, then, for — l/y/2 < 
V < l/y/2, 



P{Y < y) 



dp 



h-H[-oo,y}) 



1 



cos fc/\/l+cos 2 k<y 
1 



arccos 



7T 



dk 

47T 

k 



which gives as density f(y) of Y, 

dy 



f(y) 



7r(l-^)v/r^2? 



(35) 



in agreement with the result of |2j . The same 
holds for every unbiased walk defined by (|23|) . 

The above result can be interpreted as the 
weak convergence of the sequence X n /n of 
operators on H, as n — > oo, to an operator 
V, defined on a dense subspace of H with 
spectral resolution 



V 



E 



D\j{k) 



dEj(k), (36) 



It follows that [min h, max h] = [— Jp~, J~f\i where dEj(k) is the projector over the 
whence the bias factor of the walk sets a limit eigenspace corresponding to the eigenvalue 
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Xj(k) of U(k). (The weak convergence of un- 
bounded operators here is formally defined as 
the weak convergence of the corresponding 
unitary operators exp(isX n /n) — > exp(isV) 
for every real s.) The limit operator is diago- 
nal in the eigenbasis of the unitary evolution 
of the walk and gives 

(X n ) ~ (V)n, (37) 

that represents the Heisenberg equation of 
motion for the position, in the limit n — > oo, 
if we interpret V as the "velocity" operator. 
Thus, asymptotically, the centre of the wave 
packet moves with constant speed, given by 
V. It is worth pointing out that, although 
the equation of motion resembles the one of 
a classical system with constant velocity, the 
state of the quantum particle spreads in time, 
with a quadratic growth in the variance of the 
position distribution. 

III. WEAK LIMIT FOR 
(f-DIMENSIONAL QUANTUM WALKS 

Let d > 1. The classical random walk 
on the integer lattice Z d models the mo- 
tion of a particle that moves in an unbiased 
manner in a d- dimensional space. Let e^, 
i G {1,2, ... ,n}, be the unit vector in the 
direction of increasing zth coordinate. Let 
Ri,R.2,... be independent identically dis- 
tributed random variables, each being uni- 
form on the set {±ej : i = 1, 2, . . . , d}. The 
position of the particle at time n is given as 
the sum 

n 

X n = ^% (38) 

3=1 

By the central limit theorem for d- 
dimensional random walk, the random vec- 
tor ~X. n /yfn converges weakly as n — ► oo to 
a random vector in M. d having the multivari- 
ate normal distribution N(0, 1/d), where I is 
the d x d identity matrix. We shall see in 
the following that a corresponding weak con- 
vergence holds for a d- dimensional quantum 
random walk. 



The 1-dimensional quantum random walk 
of the last section may be extended to d 
dimensions as follows. Let ei,...,e 2d de- 
note the 2d possible shift vectors ±e$, i = 
1,2, ... ,d. The state of the system is a vec- 
tor * = {^{k)j)f =1 G H = L 2 (K d ) <g> C 2d 
where k = [k\, k 2 , . . . , kj) and the Jth com- 
ponent corresponds to a shift by the vector 
ej. At each time, the state is transformed by 
applying a rotation A acting on C 2d , followed 
by a d- dimensional shift on L 2 (K d ), cf. (HOIK 

S (c %(k)j = e i6j - k ^(k)j. (39) 

The general unitary operator that evolves the 
walk from time n = is thus 

U(k) = V{e iei - k , e ie ^}A, (40) 

where D denotes the 2c? diagonal matrix. The 
operator U (k) can be diagonalized in H, and 
has 2d eigenvalues and 2c? eigenvectors. As- 
sume that one may choose the latter as C°° 
functions of k. (See the remark at the end 
of the section.) Let fj(k), Aj(k) be re- 
spectively the eigenvectors and eigenvalues of 
U (k), with J = 1, 2, ... , 2d. The initial state 
of the system can be written in this basis as 

2d 

*o(k) = ^Mk),*o(k)Mk) (41) 
j=i 

and the state at time n as 

2d 

*„(k) = 5^A}(k)(t; J (k),*o(k)>t; J (k). 

j=i 

(42) 

The d- dimensional position opera- 
tor X® = {X 1 ,X 2 ,...,X d ) acts on 
L 2 (K. d ) as the differential vector operator 
= (—id/dki,—id/dk2,...,—id/dk d )- 
By considering each component of 
separately, it is easily seen that the operators 
X^ n converge weakly on H, as n — ► oo, to 
the corresponding components Vi, where 
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where dEj(k) denotes again the projector for all c = (ci, c 2 , . . . , q) G then X. 
onto the eigenspace of t/(k) with eigenvalue Y. 
Aj(k). This does not imply, however, that 
the sequence of random vectors associated 
with the process converges weakly on Q = 
K rf x {1, . . . , 2d}. In general, the evolution op- 
erator U (k) generates entanglement between 
the different spatial directions and it is neces- 
sary therefore to consider also the correlation 
terms between different components of X^ d \ 
The so-called Cramer-Wold device en- 
ables a simplification: in order that a se- 
quence of random variables converge weakly, 
it suffices that all linear combinations con- 
verge. More properly, we have the following, 
see p, Theorem 29.4. 



Suppose for simplicity that d = 2. For 
fixed r, we compute the expectation 



E 



d=i 



Theorem 2. Consider a sequence X n = 
(X hn , X 2 , n , ■ ■ ■ , X d>n ) , n > 1, of random d- 
vectors, and let Y = (Yi, Y 2 , . . . , Yd) be a ran- 
dom d-vector. If 



p=0 



p. 



^r—p 



(45) 



5> X 

3=1 



^^CjYj as n 



oo, 



(44) 



where we have used the fact that operators 
along different directions commute. We have 



3=1 



b&m = ;>>) p Ar p (k) (D 2 \j(k)r ( Vj (k), MVMk) + ok- 1 ) (46) 
j 

D[- p [D p 2 * n (k)] = J2(n) r XT r (k) (D 1 X J (k)) r - p (D 2 Xj(k)f (vj(k), * (k))^(k) + O^ 1 ). 

(47) 



Thus, as n — ► oo, 



E 



(c x Xi + c 2 X 2 ) n 



n 



Qcr^ik, jy- p h 2 (k, jyj |( Vj (k),* (k)>| s 

£{ Cl /n(k, J) + c 2 h 2 (k, J)} r \(vj(k), * (k))| 5 

I 



dk 



, dk 

(2^ 



where /ij(k, J) = A J (k)- 1 AA J (k), i = 1,2. With Q = K 2 x {1, 2, 3, 4}, and Z n = ciXi, n + 
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c 2 X 2 ^ n , we have 
E {(Z n /n) r ] - 



(ci/ii + c 2 h 2 ) r d/i, (48) 



where /i is the probability measure on Q given 
by 

dk 

^=|Mk),vi/ (k))| ? 



on K x {J}. 

(49) 

By the method of moments as in the one- 
dimensional case and the Cramer- Wold de- 
vice (Theorem |2J), we obtain a generalization 
of Theorem ^ to the two-dimensional case. 

As a simple example, consider the two- 
dimensional generalization of the Hadamard 
matrix given by 

/ 1 1 1 1 \ 
1-1 1-1 
1 1-1-1 
1-1-1 1 



2 



V 



(50) 



In the above notation, the unitary operator 
that evolves the walk is represented by 



U(k) = V{e 



iki 



}A (51) 



The operator U(k) may be expressed thus as 
a tensor product of two one-dimensional op- 
erators that describe Hadamard walks along 
the directions defined by k + = (ki + k 2 )/2 
and k~ = (ki — k 2 )/2: 



U(k) = U(k + ) ® U(k' 



(52) 



Its eigenvalues and eigenvectors are products 
of those of U{k + ), U{k~) respectively, and 
therefore 



A i+ (* + )Aj_(Ar), (53) 

j = 1, 2, are given by and 
J = 1, 2, 3, 4 labels the pairs (j + , j_) in some 
order. Thus 



Aj(k) 
where the A, 



hi(k, j) = Xj(ky 1 D, l x J (k) 

cos(fc + ) 



± 



2^2 - sin 2 (A;+) 
A-l cos ( k ~ 



± -1 



2^2 - sin 2 (k- 



(54) 



for i = 1,2. The limit velocity operator V is 
given by V = (Vi, V 2 ), with 

Vi= fj2 h i( k ,J)dEj(k), J = 1,2,3,4. 
J J 

(55) 

The result may be extended to arbitrary 
dimension d > 2 using the same argument, 
yielding the following result. 

Theorem 3. For the d- dimensional quantum 
random walk, 



n 



(ht(Z),...,h d (Z)), (56) 



where Z is a random element of Q = K 2 x 
{1, . . . , 2d} with distribution \x given by 
and hi{k, J) = A J (k)- 1 A.Aj(k). 

The limit observable is again diagonal in 
the eigenbasis of U(k) and represents the ve- 
locity for n — > oo. 

Technical remark. We assumed above that 
the eigenvectors of U(k) can be chosen as 
C°° functions of k. We do not know if this 
is always possible when d > 2, but it can be 
replaced by the following, weaker hypothesis, 
which we believe always holds: There exists 
an open subset O of K d with full Lebesgue 
measure (that is, the complement is a null 
set) such that the eigenvectors (and thus the 
eigenvalues) can be chosen infinitely differen- 
tiable in O. (For example, this holds if there 
is any point k where U(k) has distinct eigen- 
values, because we then can choose O as the 
subset of (0, 2n) d where the discriminant is 
non-zero; we omit the proof that this set has 
the required properties.) 

Under this assumption, the argument 
above holds for every initial value that is an 
infinitely different iable function with support 
in O. (The function h will be defined on 
O x {1, 2, ... , 2d}, but that is enough.) Such 
functions are dense in H, by a standard L 2 
result. Hence, given any initial state ^o, and 
an e > 0, we can find an initial state ^ with 
II - *nll < e for which 



-X! 



n 
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Since the evolution operators are unitary, we 
have ||\I>„-\I/ £ J| = ll^o -*q|| < e for every n, 
and it follows that for any observable event A, 
the probabilities P(X n G A) and P(X^ G A) 
differ by at most 2e. Similarly, it is easy to see 
from that \P(Z G B) - P(Z e G B)\ < 2e 
for every B C f2. It is now easy to inter- 
change the two limits e — > and n — > oo and 
obtain ()56|) ; see 0|, Theorem 4.2. Theorem|3] 
thus holds for every initial state, also under 
the weaker assumption. 



shifts by any given finite set {ei, e 2 , . . . , e^} 
of vectors in Z d . The coin flip is now repre- 
sented by a unitary matrix A in C . Theo- 
rem El extends to this case, with 2c? replaced 
by N, by the same proof. 

An interesting example is when the shift 
vectors are the 2 d vectors in {—1, l} d ; thus 
each coordinate is shifted by ±1 in each step. 
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